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I. INTRODUCTION 

The SVD-2 is universal experimental setup 
(Figure 1) [l| with the following elements: Ac- 
tive Target (AT) with microstrip silicon plates 
(40 channels of electronics) and passive plates 
(C, Pb), Microstrip Vertex Detector (MVD) 
(8.5 thousand channels), Magnetic Spectrom- 
eter (MS) with Wire Proportional Chambers 
(MWPC) (near 18 thousand channels), Thresh- 
old Cherenkov Counter (TCC) (32 photo mul- 
tipliers), Scintillation Hodoscopes (SH) (on 12 
strips in vertical and horizontal planes) and the 
Detector of Gamma quanta (DEGA) with radia- 
tors from lead glass on 1344 channels. The data 
taken run was performed in the proton beam of 
IHEP accelerator with E p =70 GeV. The total 
statistics of 5 * 10 7 inelastic events has been ob- 
tained. 

For last 20 years the obtained number of 
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events with open charm production in a proton- 
nuclear interactions in experiments is much 
less than statistics of experiments with elec- 
tron beams in which the basic properties of the 
charmed particles (mass, branching of decays, 
etc.) were studied. But the data of proton- 
nuclear experiments are important for studying 
of dynamics of charmed quarks production in col- 
lisions of nucleons and their hadronization mech- 
anisms, for checking predictions of existing theo- 
retical models. 

Taking the cross section of inelastic pp- 
interaction {(Ji n {pp) = 31.44 mb at 70 Gev [2|) 
and using the experimental fact, that the cross 
section of cc-pair production in a nucleon-nuclear 
interactions linearly depends on atomic mass 
number of target, whereas "the usual" inelastic 
section is proportional A ' 7 , it is possible to esti- 
mate the total number of events with charm in ex- 
periment E-184 for statistics of 52 million events 
with inelastic pA-interactions as 

N(cc) = N * (cx(cc) * A 1 ) / (Gi n (pp) * A - 7 ), 
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SVD-2 layout 
CI, C2 - beam scintillation and Si-hodoscopes; 
C3, C4 - Si active target and vertex Si-detector (AT and VD); 
1, 2, 3 - block of minidrift tubes (MD); 

4 - MWPC of magnetic spectrometer (MS); 

5 - threshold Cerenkov counter (CC); 

6 - scintillation hodoscope (SC); 
7 -detector of -quanta (DEGA) 



FIG. 1: SVD-2 setup layout. 

where Nq - number of inelastic events in target 
with atomic mass number A, a(cc) is accepted as 
1 fib. 

In Table I the experimental number of events 
in different plates of the target and number of 
expected events with charm are presented. 

Using yields of particles and branching of 
their decays, it is possible to estimate number 
of events with different modes of charm particles 
decays (Table II). 

The program FRITIOF 7.02 [3J for the sim- 
ulation of production processes has been used for 
D decays selection. Some parameters of model 
were adjusted with D° and D spectra measured 
in OPAL [4] and CLEO @] experiments. Other 
ones were used by default. Simulation of charm 
decays products registration in SVD-2 setup was 
made with GEANT3.21 [6]. 

II. SIMULATION OF REGISTRATION 
FOR EVENTS WITH V°. 

The position of detectors and passive el- 
ements of the setup has been simulated using 
drawings with real metrological measurements 
and with results of alignment. The experimen- 



tal map of magnetic field has been applied too 
The probability of interaction in target plate 
has been calculated with taking account of its 
thickness and nuclear length. The sum of seven 
probabilities are normalized to 1. Then the plate 
number with interaction was defined with ran- 
dom number generator [0,1]. The coordinate Z 
of the interaction along beam was defined with 
the center of the plate plus uniform displacement. 
Transverse coordinates (X, Y) of the vertex were 
pointed by the experimental beam profile. Parti- 
cle kinematics in interaction point was given by 
FRITIOF7.02. There were 3 files with carbon, sil- 
icon and lead interactions. The decays of unsta- 
ble particles were realized in GEANT. Appointed 
modes of decay were used for charmed mesons. 
During particle tracking in MVD the charge dis- 
tribution along strips, the noise in each channel 
and the experimental amplitude cuts were taken 
into account. During the hit forming in MS ex- 
perimental efficiencies of MWPC were used. 

There are the following procedures of data 
handling system: 

• Filtration of MVD data and the selection 
of events with second vertex near interac- 
tion point as the candidates for charmed 
events. The analysis in track parameters 
space {a,b} [9] was used for that; 

• Geometrical reconstruction of tracks in MS 
and defining of momentum for charged par- 
ticles; 

• The events analysis, including the kinemat- 
ics of their. 

The results of simulation were checked by means 
of comparison the distributions for MC and ex- 
perimental data. There is a good consent in num- 
ber of events for each plate of the AT. Multiplicity 
of charged particles in primary vertex (Figure 2a) 
and their momentums (Figure 2b) are in a good 
consent with model too. 

i^Tg-mesons registration. 

Simulation of K® registration in pA- 
interactions is a checking procedure for the 
searching of neutral D-mesons, because cross sec- 
tion of their production is known at our energy 
from MIRABELL data at 69 GeV/c [10]. There 
were 540000 MC events with 124000 Jf° and 
85000 of them were with (ir + ir~) decays. Next se- 
lection criteria have been applied for these events: 
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TABLE I: Number of expected events with charm at a(cc) = lfib 



The 


Thickness 


N int 
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A 0.7 


N 


N(cc) 


target 


(mkm) 


(%) 






(million ev.) 


(ev.) 
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540 


21 


12 


5.7 


10.92 


732 


Si 


300*5=1500 


55 


28 


10.3 


28.60 


2472 


Pb 


270 


24 


207 


41.8 


12.48 


1966 


TOTAL: 




100 






52.00 


5170 



TABLE II: Number of expected events with D° decays 



Decay 


Branching 
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Lead 


Ntot 
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N 
events 


Yield 


N 
events 


Yield 


N 
events 


events 


D° -> K-7T+ 


0.038 


0.488 


14 


0.497 


47 


0.527 


39 


100 


D° -> K+ir- 


0.038 


0.590 


16 


0.585 


55 


0.578 


43 


114 


TOTAL: 






30 




102 




82 


214 



• the number of charged tracks > 3; 

• the distance between primary vertex and 
Kg decay > 0.5 mm; 

• Z-coordinate of K® decay < 35 mm; 

• tracks from K® decay should cross the last 
plane of MVD. 

Only 2674 (=3%) kaons with (ir + ir~) decays are 
satisfied these conditions. 

The criteria of V° events selection have been 
optimized at the simulation. They are the follow- 
ing: 

• the track from V° should have the impact 
parameter to primary vertex, 6/cr& > 2; 

• two tracks from V° are crossing in the same 
space point, x 2 < 4 [111 ]; 



• the V0 vertex is separated from primary 
one, (Z 2 - Zx)/sqrt(af + of) > 3; 

• the V° is coming through modified 
Podolanski-Armanteros criteria [l~2[ ]. 

Fitting effective mass spectrum for (-7r + 7r - ) sys- 
tem of MC events with Lorenz function for the 
signal and 2n-polynom for background we get K® 



mass = 497.5 MeV (in PDG 497.6 ), the width of 
signal =7.4 MeV. The efficiency for events with 
K° s is 0.36%. 

D° and D° registration. 

To optimize reconstruction procedure for 
two-prong decays (Kir) of neutral D-mesons 100 
thousand MC events with charm have been used. 
Number of D° was equal 51133, that correspond 
to yield of D°. The simulation allows to esti- 
mate the efficiency of D-mesons registration (it is 
7.2%) and to optimize criteria of events selection. 
In Figure 3 it is shown as efficiency of registration 
depends on momentum and Xf of mesons. For 
MC events D° mass is 1864 MeV and width of 
signal is 33 MeV. The registration of two-prong 
decays of neutral anti D-mesons has been simu- 
lated too (Figure 4). In the aperture one gets 23% 
D° only because they fly back in c.m.s. mainly. 
Total efficiency of the registration is 2.7%. For 
MC events D° mass is 1866 MeV and width of 
signal is 36 MeV. 



III. EXPERIMENTAL DATA. 

Registration of ^-mesons. 

To compare the results of simulation with 
experimental data, 1115091 experimental events 
with interaction in the active target have been 
processed by the same procedures. The fitted sig- 
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FIG. 2: FRITIOF predictions (line) and experimental 
data (• - C,M - Si, A. - Pb): 

a) the multiplicity of charged particles in primary ver- 
tex; 

b) the momentum of charged particles. 



nal gives value of 498.6 MeV for K° s mass, width 
of signal is 12 MeV, number of events in the sig- 
nal is equal 222±20. In Figure 5 momentum (a) 
and Xf (b) of reconstructed into strip of mass 
spectrum signal after leveling of histograms on 
the number of entries are compared for MC and 




FIG. 3: The momentum (a), Xf (b) for all and recon- 
structed D° (shaded); the efficiency of D° registration 
depending on momentum (c) and Xf (d). 




FIG. 4: The momentum (a), Xf (b) for all and recon- 
structed D° (shaded); the efficiency of D° registration 
depending on momentum (c) and Xf (d). 



experimental events. 

Selection of events with feasible decays 
of D°-mesons. 

To select the events with possible two-prong 
decays (Ktt) of neutral D-mesons next basic cri- 
teria were applied: 

1. The distance between primary vertex and 
V° decay should be more 0.5 mm; 

2. The tracks from V° decay should have the 
impact parameter to primary vertex, but 
V° track should be pointed to it; 

3. The effective mass of system (Ktt) should 
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FIG. 5: The momentum (a) and Xf (b) of recon- 
structed for MC events (line) and for experimental 
events (points). 



be in the range ±0.5 GeV from D° mass 
(=1.865 GeV); 

4. The momentum of system (Kir) should be 
more 10 GeV/c; 

5. Pt of decay particle should be more 0.3 
GeV/c to moving {Kit) system according 
to Podolanski-Armanteros criteria and to 
suppression condition for background from 
kaon and A (Figure 6); 

6. The accepted hypothesis from two {K~tt + ) 
and (K + tt~) modes has the mass nearest to 
D° mass; 

7. V° should satisfy to selection criteria dur- 
ing visual inspection of events by physicist. 

It has been analyzed near 700 events. Prin- 
cipal causes for event rejection - the presence 
of background particles which were not recon- 
structed in space, but can belong to secondary 




-I -0,9 -0.6 -OA -0.2 
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FIG. 6: Podolanski-Armanteros plot for K®, A° and 
D°. 



vertex, the track from V° can belong to primary 
vertex, the possibility to accept V° instead of sec- 
ondary interaction in the next plate of the target, 
etc. Additional background sources are: the ig- 
norance of real alignment and magnetic field, the 
intrinsic background from charm decay modes. 

The effective mass spectra (Ktt) system be- 
fore (a) and after (b) checking of events by physi- 
cist are presented in Figure 7. Because of small 
statistics, spectra for (K~n + ) and {K +/ k~) sys- 
tems were united in one spectrum. The area of 
mass, for which physical examination was carried 
out, has been limited from 1.7 to 2.0 GeV. In 
peak the signal/noise ratio is (51±17) / (38±13). 
Fit by straight line plus Gauss function for data 
after checking of events by physicist gives for 
D° mass 1861 MeV and width of signal a=21 
MeV. Parameters of (Kir) systems for the se- 
lected events in Mo ± 3a range (momentum, Xf 
and decay length) for MC events (line) and ex- 
periment (points) are presented in Figure 8. It 
is visible that the measured properties of exper- 
imental (Kir) systems correspond to the proper- 
ties of D°-mesons from simulation. 



IV. ESTIMATIONS FOR CHARM 
PRODUCTION CROSS SECTION. 

To estimate the cross section, except knowl- 
edge the number of events in the signal and the 
total statistics, the knowledge of other factors is 
necessary too: trigger factor, efficiencies of all 
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FIG. 8: The decay length, the momentum and Xf of 
(Kir) system for MC and for experimental events. 



FIG. 9: Charge particles multiplicity in the primary 
vertex for MC (line) and experimental (points) inelas- 
tic pA-interactions. 



procedures of data handling system. 

The trigger factor, i. e. suppression of in- 
elastic events registration during data acquisition 
has been estimated by comparison of charged par- 
ticles multiplicity in primary vertex for MC and 
experimental inelastic events. The analysis of 
these distributions (Figure 9) gives K tr i g = 0.51. 

The fact is that the efficiency of the handling 
procedures decreases with experimental data. 
This factor is hardware factor and it can be es- 
timated using of signal. production cross 
section at our energy for pp-interactions is known 
(<7=3430 lib) and A-dependence of cross section 
for pA-interactions is A a , where a=0.78 [l~3j ]. The 
efficiency for MC events with K® is 0.36% and the 
number of decays K% — > it + it~ may be 

N(K°) = 52 * 10 6 * (3.43/31.44) * (A - 78 /A - 7 ) * 
0.692 * 0.0036 = 19200, 

as average atomic mass number of target nu- 
clei is A=68. We have experimental signal 
from K® equal to 12000 decays only. Then 
^=19200/12000=1.6. 

The estimation of charm production cross 
section in our experiment has been made by two 
ways. We have: 
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• Branching of decay D° — > Kir: Br =0.038 


• Yield of D° = 49% [3 

• Number of events in signal: 
N reg = 7(C)+22(Si)+22(Pb)=51 

• The registration efficiency (average on hy- 
pothesis): e = 0.036 

• Predictions for the number of events with 
a(cc) = l fib: 

N P red= 30(C)+82(Si)+102(Pb)=214 

• Trigger factor K tr i g = 0.51 

• Hardware factor K^d = 1-6 

• Integrated luminosity: L, ni = N p A/e/<J p N 

= [0.48(C)+1.27(Si)+0.55(Pb)] *10 33 =2.3*10 33 
Method 1. 

The first calculation of cross section is based 
on a prediction of the number of D° decays in our 
statistics from the assumption of linear depen- 
dence of charm cross section on A and ^4°' 7 for 
cross section of inelastic interactions. Though, 
number of the allocated decays is small to leave a 
problem of averaging of atomic mass number, the 
cross section was calculated separately for each 
material of a target. 

a(cc) = K hd * (N reg /e)/(N pred / K trig ) = 
= 5.3 (C); = 4.9 (Si); = 6.1 (Pb). 

And the weighted average cross section is 

a(cc)= 5.5 (fib I nucleon) . 



Method 2. 

In the second calculation 
nosity from the paper 



15] 



integrated lumi- 
is used. Integrated 
luminosity for each material of a target is de- 
fined, multiplying it by factor from Table I. As 
this luminosity corresponds to pA-events its value 
needs to be divided by A that corresponds to 
A-dependence of cross section for inelastic inter- 
actions. We are obtaining D° production cross 
sections on a nucleus for different materials: 



&nuc(D°) = K hd * N reg /(Br * e)/(L int /A - 7 ) = 
=96.5 (C); = 209.5 (Si); = 1949.0 (Pb) (fib) 

These values of cross sections show depen- 
dence on atomic mass number of target nuclei 
(A-dependence) with parameter a = 1.08 ± 0.12. 

To compare these values to the cross section 
found above (Method 1) and to the data of other 
experiments, it is necessary to take into account 
A-dependence of charm production cross section 
(a=1.0) from other experiments. Then we have 
D° production cross sections on a nucleon: 



a(D°) 



mc (D°)/A = 8.0 (C); = 7.5 (Si); 
9.4 (Pb) (fib / nucleon) 



To estimate value a(cc) it is necessary to di- 
vide this value on production yield of neutral D- 
mesons which is equal 49% by results of measure- 
ment in experiment [141 ] and on number 2 since 
we consider the sum (D° + D°). Then we have 

a(cc) = a(D°) / 0.49 /2. = [8.2 (C); = 7.6 (Si); 
= 9.6 (Pb) ] = 8.7 (ub/nucleon) 

We estimate uncertainties of measured cross sec- 
tion as 34% (statistical) and 20% (systematic). 
And the total error is 54%. 



V. CONCLUSION. 

The result of data handling in E-184 exper- 
iment for selected neutral D-mesons decays gives 
the estimation of charm production cross section 
in pA-interactions at near threshold energy 70 
GeV: 

a(cc) = 7.1 ± 2.4(stat.) ± 1.4(syst.) 
(fib /nucleon). 

The dependence of cross section versus energy in 
c.m.s from fl6l ] is presented in Figure 10, includ- 
ing this result. The theoretical limits for cross 
sections are presented too. 

The attempts for charm cross section esti- 
mation at threshold energy were made 20 years 
ago with BIS-2 IHEP setup, when carbon tar- 
get was irradiated by neutrons at 40-70 GeV 
[171 ] . In kinematical range Xf > 0.5 measured 
D° production cross section occurred much larger 
of theoretical predictions and is o~(D°) = 28 ± 
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14(ab/ nucleus). The cross section in our exper- 
iment is larger then theoretical estimations too. 
This fact demands more detailed study and much 
bigger stastics. The executed work is important, 
because the information on charm cross section at 
near threshold energy in pA-interactions is poor 
enough now. 
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FIG. 10: Experimental cross section of D° production 
in pA-interactions with the results of E-184 experi- 
ment. 
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